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Synchronization vs. decoherence of neutrino oscillations at intermediate densities 
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We study collective oscillations of a two-flavor neutrino system with arbitrary but flxed density. 
In the vacuum limit, modes with different energies quickly de- phase (kinematical decoherence), 
whereas in the limit of infinite density they lock to each other (synchronization). For intermediate 
densities, we find different classes of solutions. There is always a phase transition in the sense of 
partial synchronization occurring only above a density threshold. For small mixing angles, partial 
or complete decoherence can be induced by a parametric resonance, introducing a new time scale to 
the problem, the final outcome depending on the spectrum and mixing angle. We derive an analytic 
relation that allows us to calculate the late-time degree of coherence based on the spectrum alone. 

PACS numbers: 14.60.Pq, 97.60.Bw 



I. INTRODUCTION 

Flavor oscillations in dense neutrino gases exhibit col- 
lective phenomena caused by neutrino-neutrino interac- 
tions (ll-[l7t. Recently this topic has been studied in- 
tensely in the context of supernova neutrino oscillations; 
for a recent review see Ref. [l^. In much of this liter- 
ature, numerical tools were used to study cases of prac- 
tical interest for supernova neutrinos. Moreover, there 
has been progress in the analytic understanding of such 
intriguing phenomena as spectral splits, caused by the 
adiabatic decrease of the neutrino density with radius in 
a supernova. However, our theoretical understanding of 
many aspects of collective oscillations remains unsatis- 
factory. This applies, in particular, to the question of 
angular synchronization vs. angular decoherence of neu- 
trinos streaming from a supernova core [TT| - [l3j . 

With such general questions in mind, we here return 
to an elementary case of collective effects: synchronized 
oscillations in an isotropic neutrino gas with fixed den- 
sity. Modes with different energies, E, and thus differ- 
ent vacuum oscillation frequencies, w — Am^/2i? with 
Am^ the neutrino mass-squared difference, quickly de- 
phase (kinematical decoherence). However, if the density 
is large one obtains synchronized oscillations: all modes 
"stick together" and oscillate with a common frequency 
Wsync = (w) 0- In the language of flavor-polarization 
vectors P;^, all modes together form a common global 
P = / duj'Puj that precesses as a single object. 

While the extreme cases of very dilute and very dense 
gases are easy to understand, the question of what ex- 
actly happens for intermediate densities is less obvious. 
Different possible answers come to mind, (i) The system 
could decohere completely, but perhaps on a time scale 
that grows longer for larger neutrino densities, (ii) A sta- 
tionary state with some global P of reduced length could 
be reached (partial decoherence). 

The correct answer has elements of both ideas, in detail 
depending on the neutrino spectrum, density, and mix- 
ing angle. In the simplest case of maximal mixing, the 



system decoheres completely for densities below a certain 
threshold value. For larger densities, synchronization is 
partial and becomes perfect in the limit of infinite den- 
sity. The behavior of the system as a function of density 
thus shows a behavior similar to a phase transition. 

Other forms of behavior arise for smaller mixing angles. 
The main reason is that polarization vectors that are ini- 
tially aligned with each other and almost aligned with 
the effective magnetic field responsible for their preces- 
sion can later deviate strongly from both directions due 
to parametric resonance effects. Depending on whether 
or not the resonance frequency falls within the spectrum, 
complete decoherence can arise, but can take a very long 
time to reach equilibrium. In other words, the decoher- 
ence time scale can be very different from any other time 
scale of the problem. 

Our study is in many ways complementary to Panta- 
leone's investigation of the Stability of incoherence in an 
isotropic gas of oscillating neutrinos 0]. He asked if a 
completely incoherent initial state would spontaneously 
develop global flavor polarization, a process resembling 
the spontaneous synchronization of a collection of cou- 
pled oscillators. He found neutral stability, implying that 
spontaneous synchronization does not occur, largely be- 
cause the equations of motions conserve energy. This 
contrasts with some models designed to mimic sponta- 
neous synchronization of certain biological systems pj . 

Returning to the polarization decay of an initially syn- 
chronized neutrino ensemble, we begin in Sec. H] by set- 
ting up the equations of motion of a system of interacting 
flavor polarization vectors. In Sec. lIIII we consider simple 
cases of kinematical decoherence. In Sec. lIVI we derive an- 
alytic relations that allow us to calculate the asymptotic 
behavior. In Sec.|V]we show that the threshold behavior 
of partial synchronization can be understood if one rep- 
resents the spectrum in terms of two discrete polarization 
vectors. In Sec. IVII we identify the phenomenon of self- 
induced resonant decoherence. In Sec. lVIll we summarize 
our flndings. 
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II. EQUATIONS OF MOTION 

A. Matrices of occupation numbers 

For the kinetic evolution of a particle ensemble, the 
fundamental entities are occupation numbers /p = 



(a^ftp) for every mode p. Here aj,ap is a number op- 



erator and (■ ■ • ) an expectation value or ensemble aver- 
age. Including the flavor degree of freedom, the occu- 
pation numbers are replaced by flavor matrices of occu- 
pation numbers gp,ai3 = (Op.aOp,/?) where a and /3 are 
flavor indices. The diagonal entries represent the occu- 
pation numbers for different flavors and are related to 
the survival probabilities, whereas the off-diagonal ele- 
ments encode phase information [l^, [2d| . Based on this 
formalism, one can derive a Boltzmann collision equation 
that allows one to study the evolution of a neutrino en- 
semble under the influence of both flavor oscillations and 
collisions 

Here we study a much simpler case where collisions are 
neglected, motivated by neutrinos streaming freely from 
a supernova core. However, in our schematic model we 
ignore all geometric effects and consider the time evolu- 
tion of an isotropic ensemble rather than the spatial vari- 
ation of a stationary flux. In the limit of ultrarelativistic 
neutrinos, the equations of motion (EoMs) are [1^ 



AP r- 

2^ + ^^^ 



F g, Qp 



(1) 



where [ • , • ] is a commutator and Af ^ is the neutrino 
mass-squared matrix. It is not diagonal in the weak in- 
teraction basis, causing vacuum flavor oscillations. The 
global matrix ^) = / gp d'^p/(27r)'^ has diagonal entries of 
neutrino number densities. This "neutrino matter term" 
has the same structure as the ordinary matter term where 
the role of the global g is played by the matrix of electron, 
muon and tau-lepton number densities that is diagonal in 
the weak interaction basis. Indeed, a background of neu- 
trinos would be fully analogous to ordinary matter, were 
it not for the off-diagonal elements of g as first recognized 
in Pantaleone's seminal paper 

Full kinematical decoherence corresponds to the off- 
diagonal elements of g vanishing in the mass basis. At 
a sufficient distance from a source, a neutrino flux will 
always end up in such a state. It is sometimes described 
as the wave packets of different mass eigenstates no longer 
overlapping. 



B. Polarization vectors 



and 8-dimensional polarization vectors). Likewise, we 
may write M^/2E = ^Tr{AP/2E) + iwB • cr where 
u! — Am?/2E is the vacuum oscillation frequency and 
B a unit vector, playing the role of an effective magnetic 
fleld in the spin-precession formula below. 

We study an isotropic ensemble, i.e. we consider modes 
gE that represent all momentum modes of the same en- 
ergy i5 = |p|. It proves more convenient to label the 
modes with cu — Am^/2E and corresponding polariza- 
tion vectors P^. The EoMs then appear in the familiar 
precession form (Fig. [1]) 



Huj X Puj where Huj — wB + /iP . 



(2) 



The total polarization vector is P = / duj P^ which we 
normalize as |P| = 1. In this case iJ, — ^/^Gpn^, with 
the neutrino number density, is the interaction energy 
felt by a neutrino in the mean field of all others. 

The precession of individual modes is around B which 
thus marks a natural symmetry direction of the problem. 
Therefore, we choose the B-direction as the z-direction 
of a cartesian coordinate system in flavor space. We as- 
sume that all neutrinos are initially prepared in the same 
flavor, so all P^^ initially point in the weak-interaction 
direction. It is tilted relative to B by twice the mixing 
angle 9 (Fig.[T]). Therefore, we have 



B 



and P^(0) 




(3) 



where g,^ — g{uj) — |Pi^(0)| is what we call the neutrino 
spectrum normalized to J du) gi^ — I. The initial di- 
rection in the x-y-plane is arbitrary, owing to a global 
SU(2) symmetry of the problem: the absolute phase of 

relative to another flavor i^j. is arbitrary. 

The EoMs conserve the quantity J dw wB • P^^ + ^ fiP^ 
which thus plays the role of total energy. This is seen by 




For many studies it proves more transparent to rep- 
resent gp in terms of polarization vectors Pp in "flavor 
space." We restrict our investigation to two flavors where 
gp is a Hermitean 2x2 matrix and thus can be expressed 
in terms of Pauli matrices as gp = ^TT{gp) + ^Pp ■ cr 
(while for three flavors, we need Gell-Mann matrices 



FIG. 1: Precession of a polarization vector P,^ with frequency 
uj around the z-direction (the mass direction in flavor space) . 
The initial orientation is in the weak-interaction direction, 
tilted by twice the mixing angle relative to the z-direction. 
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taking the time derivative and inserting the EoMs on the 
r.h.s. for P;^. In other words, the EoMs can be derived 
from a classical Hamiltonian. This is seen most easily 
if we use discrete polarization vectors for a set of N 
frequency bins uji where P = P^, 

N N 

i/ = ^^.B.P, + | ^ P,.P,. (4) 

i=l i-.j — ^ 

The polarization vectors play the role of classical angu- 
lar momenta with Poisson brackets [Px,Py\ = Pz and 
cyclic permutations.^ They inherit this property from 
the SU(2) algebra of the "flavor isospin matrices" [lo| . 
The precession equations follow from the classical Hamil- 
ton equations of motion P^ — [Pi, H]. The Pi are classi- 
cal variables because they arise from occupation number 
matrices which themselves are classical objects, not quan- 
tum operators, as they represent an expectation value. 

The quantity P plays the role of the total angular mo- 
mentum of the ensemble. Naturally, its component along 
B is conserved. The conservation of B • P is also easily 
seen by taking the time derivative of this bilinear and 
inserting the EoMs for P^. 

In the polarization- vector language, our ensemble of in- 
teracting neutrinos is equivalent to an ensemble of classi- 
cal angular momenta P^ that couple to an external mag- 
netic field with strength uji and to each other with a 
dipole-dipole interaction of equal strength fi. However, 
they do not couple by some nearest-neighbor interaction, 
but rather every mode couples to every other with the 
same strength. 



C. Antineutrinos 

Occupation number matrices for antineutrinos can be 
defined in three different ways, all of which have been 
used in the literature. If 6j, and bp are creation and an- 
nihilation operators for antineutrinos in mode p we may 
define Qp^ap = {bl^abp^fs) or {bljsbp^a) or 1 - (bl^^bp^a), 

the latter arising from normal ordering of (&p,a^p^)- 
The third definition amounts to the flavor-isospin inter- 
pretation of polarization vectors [l^ and is the most 
useful in our context. Antiparticles are interpreted 
as negative-energy states and are thus described by 
(jj = Am^ /2E < 0. In the polarization- vector language, 
antineutrinos are then described by with negative uj. 

The fiavor interpretation is now as follows: P^j in 
the positive weak- interaction direction means i^e, in the 
negative weak-interaction direction means Vx if w > 0, 
whereas for uj < these cases mean and Dg, respec- 
tively. This sounds complicated, but is simple in the hole 



interpretation of antiparticles where the presence of in 
mode p is equivalent to the absence of a negative-energy 
j/fi from a filled mode in the Dirac sea. 

Once we are studying the dynamics of our abstract 
Hamiltonian problem of interacting spins, the detailed 
flavor interpretation is irrelevant. We simply consider an 
ensemble of P,^ with a range of positive and negative u. 

D. Rotating frames 

The dynamics of the system is the same if we go to a 
frame rotating with loq around the z-direction This 
is seen most easily if we consider the spin-spin Hamilto- 
nian of Eq. (jH) and recall that B • P is conserved. There- 
fore, we may add a term —loq'B ■ P, leading to 

N N 

iJ = ^(w.-^o)B.P, + | ^ P,.P,. (5) 

i=l i-.j — ^ 

Therefore, even if we are having in mind an ensemble 
of neutrinos (no antineutrinos) and therefore a range of 
positive w, we may shift the spectrum g{uj) by a conve- 
nient amount —ujq and thus include negative w. In other 
words, the negative-w part of the spectrum can be inter- 
preted as antineutrino modes or as neutrino modes after 
such a shift. The possibility of this general interpretation 
and the seamless connection of neutrino and antineutrino 
modes is the main advantage of the flavor-isospin inter- 
pretation of polarization vectors. 

E. Spectra 

We always consider non-negative spectra, > 0, 
avoiding "spectral crossings." The latter lead to unsta- 
ble forms of motion ( "flavor pendulum" ) that are at 
the origin of spectral splits [13] . Addressing the question 
of decoherence in such more general cases is beyond our 
present ambition. 

We will explicitly study two cases of simple spectra^ 
that are symmetric about w = and normalized to unity. 
One is a box spectrum where — everywhere except 

5,.^ = ^ for - 1 < w < +1 . (6) 
The other example is a Gaussian with unit variance 

We show these spectra in Fig. [51 



^ To simplify notation we take frequencies and time to be dimen- 
Here [ • , ■ ] is a Poisson bracket, not a commutator. Of course, sionless, or rather, to be expressed in terms of some fiducial fre- 

'yPx- quency ijJq so that ui ui/uiq, fi —> fi/uiO) and t — > coot. 
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FIG. 2: Box spectrum (red, dashed) and Gaussian spectrum 
(blue, solid) as defined in Eqs. © and © respectively. 

The box spectrum is representative for "compact" 
spectra that are non-zero in a finite w-range, whereas the 
Gaussian is representative for spectra with long tails. We 
will see that spectral tails are crucial for the phenomenon 
of resonant decoherence to be discussed in Sec. IVll 




Time 

FIG. 3: Decay of the order parameter R{t), defined in Eq. ((8|, 
for the box spectrum (red, dashed) and the Gaussian spec- 
trum (blue, solid) as given in Eq. (llOp . 

neighboring modes, leading to a shrinking of after 
coarse-graining over small but finite frequency bins Aw. 
For our simple cases we are using a single frequency bin: 
the entire spectrum. 



III. KINEMATICAL DECOHERENCE 



B. Decoherence in the non-interacting case 



A. Order Parameter 

When neutrino modes de-phase relative to each other, 
the system becomes in some sense disordered, but it is 
not trivial to define an objective measure for this effect. 
Flavor oscillations alone (no collisions) do not increase 
the entropy of the system Kinematical decoherence 
is reversible. For example, a neutrino flux decoheres kine- 
matically after a short distance from the source, but a de- 
tector with sufficient energy resolution can still pick out 
oscillations in an energy-dependent way and in this sense 
restores coherence. Another example is the phenomenon 
of spin echo: reversing the direction of the magnetic field 
unwinds the phases accrued by different modes and the 
original polarization re-appears. 

In our case the simplest quantity to measure coherence 
is what Pantaleone calls the "order parameter" Q 



R{t) 



\P±{t)\ 
sin 20 



(8) 



where it was assumed that |P(0)| ~ 1. Here is that 
part of P transverse to B. Since B • P is conserved, 
it is only the transverse component that can shrink by 
kinematical decoherence. Dividing by sin 29 normalizes 
the initial order parameter to i?(0) = 1. 

Were we to study decoherence for more complicated 
spectra, this definition would be less useful. When spec- 
tral crossings and pendular motions are involved, Pj_ can 
grow even though the system becomes "more disordered" 
as it evolves. Probably one would need to define kinemat- 
ical decoherence in a differential sense, the dephasing of 



De-phasing is easy to understand in the dilute limit 
where neutrino-neutrino interactions play no role (^ = 0) 
and the evolution is determined by vacuum oscillations 
alone. The transverse components evolve as Px{t) = 
sm29 J dojgu, cos(a;i) and Py{t) — sin 26* Jdcjg^ sin(a;t). 
In other words, the order parameter evolves as 



m 



dcj e" 



(9) 



Therefore, R{t) is given by the Fourier transform of g^i- 
Unless guj includes "spikes" of the form S{uj—uji) at one or 
more frequencies uji, the late-time evolution is R{t) — > 
for i — )■ oo (complete kinematical decoherence). 

For our two explicit examples, the time evolution of 
the order parameter is (Fig. [S]) 



Box: R{t) 
Gauss: R{t) 



sin(i) /t , 



exp( 



(10) 



The box spectrum has "hard edges" and thus shows the 
usual "ringing" of its Fourier transform. 



C. Self- maintained coherence 

For non-zero interaction strength, the system can show 
"self-maintained coherence" in that an oscillatory motion 
persists at late times. We illustrate this in Fig. |4] where 
we show the evolution of the order parameter R{t) for 
the Gaussian spectrum with maximal mixing and dif- 
ferent values of fi. For /i above a threshold value 
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FIG. 4: Order parameter for the Gaussian spectrum with 
maximal mixing and = 0, /xo, 1, 2 and 4 (bottom to top), 
where fio ~ 0.7979 is the threshold value for partial coherence. 
Dashed lines: Solution when P(t) in the EoMs is replaced with 
the asymptotic solution A(t) such as if P{t) would take on its 
asymptotic solution instantaneously. 



to be identified later, R{t) at first shrinks like in the 
non-interacting case, but then oscillates with decreas- 
ing amplitude around a non- vanishing asymptotic value 
Ra- In other words, the global polarization vector P(t) 
shrinks in length and eventually reaches a nonvanishing 
asymptotic form A{t) that keeps precessing around the 
B direction with a certain asymptotic frequency uja- Of 
course, in the extreme case /z — >■ cjo we have perfectly 
synchronized oscillations and the asymptotic length is 
\A{t)\ = |P(0)| = 1 and iJA = wsync = (w). 

An asymptotic solution P{t) — )• A(t) performing a 
pure precession around B is the simplest case of collective 
motion. We may call it single-mode coherence in that the 
asymptotic solution is described by a single asymptotic 
polarization vector. 

However, we may construct more complicated spectra 
g{uj) where the asymptotic solution contains two or more 
modes. For example, if g{uj) consists of two separate box 
spectra, then for vanishing fi the spectrum completely 
decoheres and for infinite /x the system acts as one syn- 
chronized mode. However, one can choose fi so large that 
each individual box acts as one collective mode, being in- 
ternally synchronized, yet so small that each box essen- 
tially precesses freely around B. We now have effectively 
two weakly coupled modes and thus bimodal coherence. 
In the same way we can construct a hair-comb spectrum 
consisting of N boxes, leading for intermediate /i to A^- 
modal coherence. 

A priori it is not obvious if a given spectrum, for a 
given interaction strength /i, will decohere or not, and in 
the latter case, how many modes are required to describe 
the asymptotic behavior. In the rest of the paper we 
study only cases of single-mode coherence where the final 
behavior is pure precession. This is the case for both the 
box and the Gaussian spectrum. 



D. Sudden approximation 

We can qualitatively understand the behavior in Fig.H] 
if we imagine that P{t) shrinks instantaneously to its 
asymptotic form A(t), an assumption that we call "sud- 
den approximation." We can use the numerical asymp- 
totic solution and evolve the initial state with A{t) in- 
stead of P(i) in the EoMs, in this way reproducing what 
we would get if the sudden approximation were indeed 
correct. We show examples as thin dashed lines in Fig.|4l 

In the sudden approximation, kinematical decoher- 
ence is once more a passive de-phasing effect except 
that the individual modes no longer precess around 
a common B, but each precesses around its individual 
Hij(t) — uj'B + ^A{t). There is no feedback of P^^ on 
and we can determine R(t) by projecting the evolution of 
the individual modes on the x-y-plane. For a symmetric 
spectrum with maximal mixing we find 



R{t) = J dujg. 



+ LO"^ COS (Vk^ + t) 



(11) 



where k = /i A±^ which for maximal mixing is k = ^ A. 
For K = fiA — this reduces to Eq. (|9]). For small t 
we may expand the cosine in Eq. pT|) and find R{t) — 
1 — i {uj'^)t^ . The initial evolution depends only on the 
variance of the spectrum, but not on /j,. 



IV. PARTIAL SYNCHRONIZATION 

A. Self-consistency conditions 

We next study quantitatively the conditions for the 
existence of a nonvanishing asymptotic solution P(t) — > 
A(t) which we assume to be a pure precession with fre- 
quency UJA- Conditions for the existence of a nontrivial 
solution for fi < oo can be formulated if we study the 
system from the perspective of a frame that co-rotates 
with A in analogy to the approach taken to study spec- 
tral splits [m, [l6] . The assumed final steady-state EoM 
is given by Eq. ([2]) with 



11^ = {uj- uja) B + ^A , 



(12) 



where in the co-rotating frame A is a static vector. The 
evolution of every individual mode is now a simple pre- 
cession around this fixed vector and we refer to these 
asymptotic modes as A^^ (t) . 

At late times, these modes must have de-phased, mean- 
ing that their oscillation phase relative to their carries 
no correlation with each other. Therefore, the total po- 
larization vector is given by the time averages of A^ as 
A = J duj (Alj). The time-average of A^ is the projec- 
tion on its Hamiltonian vector 



H2 



(13) 
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implying 



/ 



do; 



[{lu - uja)^ + ^A] 



(14) 

Recall that the A;^ in this equation are time dependent, 
precessing around the asymptotic H^. 

The quantity P|| = cos 29 is conserved and therefore 
identical with . It will be useful to introduce 



K = fJ-A± and lu^ = uja — fJ- cos 20 , 



(15) 



where uj^ plays the role of a resonance frequency. The 
components of the self-consistency relation perpendicular 
and parallel to B become 



All = 



dU! 7 rr: ^ (U — CJr ) , (lo) 



where k and uir are the quantities to be determined. 

What is the meaning of the resonance? A given P^^ 
precesses in the co-rotating frame with the frequency 
io — UJA + M cos 29. When this quantity vanishes, P^j 
does not precess around the z-direction and only feels 
the transverse field. This is similar to spin resonance in a 
longitudinal B field and a rotating transverse field. If the 
transverse field is on resonance, the spin can completely 
flip even if the transverse field is weak. We will see that 
this resonance effect is indeed important for small mixing 
angles if Wr falls within the range of the spectrum. 



B. Sudden approximation 

The final angle between A,^ and the corresponding 
asymptotic H;^ is not known, so we need to make ad- 
ditional assumptions. Taking the non-interacting case as 
a guide, decoherence happens on a timescale of the same 
order as the individual oscillation frequencies and thus 
is not adiabatic. Therefore, we use the sudden approxi- 
mation where the initial P shrinks instantaneously to the 
final A. We will see that this heuristic procedure leads to 
an apparently very general result even though the sudden 
approximation can be badly violated. 

We thus evaluate the integrand in Eq. (IT51) at time 
t = 0, assuming that at this time A^ coincides with the 
initial Pi^ so that A^± = sin 20 and A^\\ = g^cos29. 
With the notation s = sin 20 and c — cos 20 we find 



— — / dcj ( 



(lo — CJr) C + K S 
{lo — CJi.)2 + ' 



(uj — UJi-) C + K S 
'^^ 9'^ 1 2" ^ - '^r) • 17 



One can easily show that in the limit /i — > oo the equa- 



tions are solved by Aj^ = s and loa 
the usual synchronized oscillations. 



^sync 



One can form linear combinations of these equations 
and also use 1 = / dw . One interesting form is 

^ f, {l + t^)RA f, il + t^)i 



where ^ = {uj—lo^) / ^ic = {uj—loa) //ic+1 and t = tan 20 = 
s/c. We recall that Ra — A±/ sin 20 is the asymptotic 
value of the order parameter with < Ra < 1. 



C. Analytic properties 

Perhaps the most illuminating form of the relations is 
dujg^ 



c 



dujgu 



{uj — Wr)^ + 



(19) 



(w — Wr)^ I - 

We may combine them to a single equation by adding the 
first, multiplied with the imaginary unit i, to the second. 
Using (w — cj]. + i'^) ("^ — — i'*) for the denominator it 
yields 

^ I du , (20) 



where is real and k is real and positive. 

Going one step further, we add an imaginary part to 
g{oj) and extend it to the complex plane by"^ 



G(f2) 



1 



+ 00 



diu' 



(21) 



This complex- valued spectrum in the complex plane ful- 
fills the boundary condition g{uj) = ReG(a;) and obeys 
Kramers-Kronig relations. The self-consistency relations 
are now expressed in the compact form 



A* 



-1 gi2e 



i7rG(aJi. + Ik) . 



(22) 



This analytic relation suggests that our heuristic argu- 
ment has provided a more general result than justified 
by its derivation. The predicted asymptotic solution ac- 
tually agrees with all numerical examples even when the 
sudden approximation is strongly violated. 

For our main examples, the extended spectra are found 
by direct integration 

i , /Q + r 

TT log 



Box: G{Q) = 
Gauss: G(Q) = 



2tt 
1 



/2Tr 



cxp 



1 
2 



(23) 



erfc 



where erfc(x) = 1 — erf(a;) and erf is the error function. 
The box case can be inverted, providing explicitly 

UJr +1K = -coth ( \ . (24) 



We use capital sans-serif letters to denote complex quantities. 
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D. Threshold for partial synchronization 

Even without explicit evaluation of G(f2), we recog- 
nize that a minimal interaction strength fiQ is required 
to achieve partial synchronization. We consider vanish- 
ing asymptotic A± and thus the k = limit. The second 
integral in Eq. (ITS)) becomes 



- = / dw - 



9u. 



■ Wo 



(25) 



where is the resonance frequency corresponding to 
At = and the integral is to be evaluated in the principal- 
value sense. Thus one can determine wq and use it to 
evaluate the first integral where for k — > the second 
factor becomes -k 5(uj — ujq). Therefore we find that 



sin 261 



(26) 



is required to obtain partial synchronization. For fi < 
one expects complete decoherence, although the equa- 
tions may have formal solutions with negative Aj_. 

For the box spectrum one finds an explicit result for 
any mixing angle, whereas for the Gaussian spectrum 
only for maximal mixing, 



Box: fiQ 
Gauss: 



(2/^)s, 

V27^ for s = 1 



(27) 




0.0 0.2 0.4 0.6 0.8 1.0 




0.4 0.6 

sm{2e) 



FIG. 5: Threshold interaction strength /io and threshold reso- 
nance frequency luq for the Gaussian (sohd blue) and the box 
spectrum (red dashed). 



For general s, we use Eq. ([22]) with k = and threshold 
resonance frequency to obtain the parametric solution 
so(wo) and /Lto(wo) with 



So = sin{arg[i G(a;o)]} 
kG(c^o)r^ 



Mo 



(28) 



Using the explicit G{uj) for the Gaussian it is straightfor- 
ward to plot Ho as a function of sin 29 in Fig. [SJ For small 
sin 20, the Gaussian case behaves completely differently 
due to resonant decoherence as we will explain in Sec. lVIl 
These results for < sin 26* < 1 (0 < 6^ < 7r/4) cover 
only the inverted hierarchy. In this sense it would have 
been better to consider the range < 6' < 7r/2 and use 
— 1 < cos 20 < as a mixing variable. For a symmetric 
spectrum, however, the results are the same except that 
Ur — —u)r in the normal hierarchy. 



V. ORIGIN OF THRESHOLD 
A. Maximal Mixing and symmetric spectrum 

The component of P parallel to B is conserved and 
only the transverse component decreases by decoherence. 
Therefore, to understand the origin of the interaction 
threshold for partial synchronization we study maximal 
mixing where P|| = cos 26* = and sin 20 = 1. To sim- 
plify matters further, we assume a symmetric spectrum 
f/(— w) = g{+uj). As a consequence, uja = ujr = ior 
any /i, i.e. the global polarization vector remains static 
except for its length P{t) and the problem reduces to 
finding this function. 

With UJA — and k = fiA± (which here is k = fJ-A), 
we need only one real-valued self-consistency relation 



1 



(29) 



In the limit /i — > oo we expand the second factor in 
Eq. dm in powers of a;^//z^ and find i?A = 1 - (w^)///^. 

For the Gaussian spectrum the threshold is at /io ~ 
\f2pK K. 0.7979. The spectrum has unit variance, so at 
large we have Ra = 1 — /i~^. Integrating Eq. ([^ using 
Eqs. ([22]) and ([231) with = 7r/4 and Q = in we find 



(30) 



1 = ^2 ^^Kyj^^^K^i 



We thus have a parametric solution for Ra (/i) in the form 
//(k) and RAifi) = i^/ shown in Fig. [51 For the box 
spectrum, the threshold is at /io = 2/7r w 0.6366 and 



Ra — ^ cotan(/i ^) . 
We show this result as a red, dashed line in Fig. [S] 



(31) 
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FIG. 6: Asymptotic order parameter Ra for the Gaussian 
spectrum (blue, solid) and the box spectrum (red, dashed). 



B. Discrete polarization vectors 



Q. 
> 



0.2 



0.0 



-0.2 



-0.4 



-1.0 -0.5 0.0 0.5 1.0 



FIG. 7: Potential V{P) defined in Eq. (l33ll for u = 0. |. 1. 
and I (bottom to top). 



When the system develops a collective mode, it set- 
tles only slowly to the final state of a pure precession, 
maintaining a long-lasting collective oscillation around 
the asymptotic solution. This observation suggests that 
the phase transition to a non-vanishing average (P) may 
already occur in a simpler system. To investigate this 
question we represent the spectrum schematically by two 
discrete polarization vectors, one for the negative and 
one for the positive frequency modes, respectively. To be 
specific, we use the normalized symmetric spectrum 



S{lu- 



(32) 



and maximal mixing. This coarse representation of the 
box spectrum does not allow for decoherence, yet we may 
study its evolution and ask for its eigenmodes as a func- 
tion of the interaction strength. 

The symmetries of our problem imply that P(t) has 
only one nonvanishing component P{t). From the full 
EoMs we find an equation of the form "kinetic plus po- 
tential energy is constant" ^P'^ + ViP) =0 with 



V{P) 



- p^Y - {I - p^ 



(33) 



We show this potential in Fig. [7] for several values of 

The motion starts at P = 1 with total energy 0. For 
< /i < 1 the motion oscillates between — 1 < P < -f-l 
and (P) = 0. On the other hand, for /i > 1 the sys- 
tem is trapped in the region of positive P. It oscil- 
lates in the range Pmin < P < 1 around the aver- 
age (P) = i(l + Pmin)- The turning point is Pmin = 
''^ so that 



(P) = 



(34) 



For /i < we have (P) = and thus we find the overall 
behavior shown in Fig. [5] 




FIG. 8: (P) given by Eq. (|34[) for the case of two discrete 
polarization vectors. 



The simple system of two discrete polarization vectors 
nicely reproduces the key features observed for a con- 
tinuous spectrum. In particular, the threshold behavior 
is primarily related to the global properties of the self- 
interacting system. Decoherence toward the asymptotic 
P{t) is a secondary effect. 



VI. RESONANT DECOHERENCE 



A. Vanishing mixing angle 



For non-maximal mixing, new phenomena can arise. 
We continue to use a symmetric spectrum g^j, but the 
co-rotation frequency for /i < (X) is no longer identical 

— 0. As an 



sync 



with the synchronization frequency u. 
extreme case we consider the limit of vanishing mixing 
angle. One may begin with the self-consistency condi- 
tions in the form of Eq. p^ . In the limit 6' 0, up 



to linear order in 6, we have t — tan (26') and 
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FIG. 9: Asymptotic values for global precession frequency 
UJA (blue) and order parameter Ra (red), assuming our box 
spectrum with vanishing mixing angle. The green (dashed) 
curve is Ra for majcimal mixing. 



FIG. 10: Evolution of the order parameter for the box spec- 
trum with wings of Eq. ((STj), using 6 = 0, 10"^, 3 x 10~^, and 
10^^ (top to bottom). 



c = cos 20 — > 1 and the conditions become 

1= [ dujg^ and 1= [ duj g^, — - — . 

J (W - J LO-LO^ 

(35) 

The second integral does not depend on Ra and thus can 
be solved and inverted to obtain Wr(M)- The result can 
be inserted in the first integral to find i?^(/x). For the 
box spectrum we find explicitly 

coih^ lii^^) — 1 
Wr(Ai) = — coth(^ ) and RaW = ^ • 

(36) 

These results also follow from Eq. in the s 
limit. We show ~uja{i^) and Ra{i-^) in Fig. [3] Note that 
— 1 < UJA < and in particular uja — > i^sync — for 
/i — ^ oo. For small /z the function i?^(/i) is very small, 
but there is no threshold. 

B. Self-induced decoherence 



which is essentially fixed by the main part of the box, 
falls within the wings. On the other hand, if 6 ^ 1 the 
evolution of P{t) can not at first differ from the case 
& = 0, so the impact of a small b can only manifest itself 
at late times. 

We study numerically an example where the resonance 
frequency is fixed to fall into the center of the negative- 
frequency wing, i.e. uj^ = —1.5, implying with Eq. (|36|) 
H = l/acoth(3/2) = 2/log(5) « 1.243 for 5 < 1. The 
evolution of R{t) is shown in Fig. [TOl For nonvanishing 
b it is exponentially damped, the rate depending on b. 

We also consider the late-time orientation of rela- 
tive to B. When Pj^ has shrunk to zero, the remaining 
motion of P^^ is a precession around the z-direction with 
a fixed angle. In Fig. [TT]we show a^j = Puj,±/{gu: sin26'), 
a quantity that is initially unity when all polarization 
vectors begin in the weak-interaction direction. We see 
that indeed the polarization vectors near cji- = — 1.5 have 
moved far away, whereas those with larger frequencies 
have actually been pulled closer to B. We show the re- 



Inspecting Eq. ([55]) we notice that the first integral 
diverges if > for lo — uj,-, whereas the second remains 
meaningful in the principal-value sense. In other words, 
Wr can be determined from the second equation, whereas 
the first implies Ra = 0, i.e. complete decoherence for 
any /z. This conclusion applies to any spectrum with 
infinite tails, whereas the resonance frequency is outside 
of the box spectrum for any /i. 

To understand better the impact of the resonance, we 
return to the symmetric normalized box spectrum, aug- 
menting it with small "wings" 

{(l-6)/2 for < 1 
b/2 for 1< < 2 (37) 

otherwise 

where < 6 <C 1. The self-consistency conditions imply 
qualitatively different solutions for any nonzero b if Wr, 



3 

4i 




Frequency cj 



FIG. 11: Final deviation of polarization vectors from the B 
direction in units of sin 29. The box spectrum with wings was 
used with b — 10~^. In the sudden approximation one would 
expect Oui = 1 (dashed line). 
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suit for b = 10 For other values the curve looks the 
same, except that the resonant peak is taller for smaller 
b. For other spectra (e.g. Gaussian), the curve is different 
in detail, but there is always the resonance peak. 

C. Intermediate mixing angle 

When the mixing angle is not vanishingly small, the 
polarization vectors cannot be driven arbitrarily far away 
from the B direction. For example, exactly on res- 
onance will completely reverse and then come back and 
so forth. For example, the resonance peak in Fig. [TT] 
reaches roughly up to Puj,± — ^Og^j sin 20. Therefore, if 
sin26' > 0.02, it is no longer geometrically possible to get 
that far and the resonance can not fully develop. 

The threshold condition for partial synchronization of 
Eq. (|26|) here reads sq = 7r/x6/2, assuming fi is such that 
the resonance falls into the wing. Specifically we consider 
H = 2/ log(5) sa 1.243 and b = 10"^ With these numbers 
we expect complete dccoherence for sin 20 < 0.01952, in 
perfect agreement with our simple geometric estimate. 

The different forms of evolution are illustrated in 
Fig. [121 showing the early evolution for the box with 
wings {b = 0.01 and sin 26 = 0.04) in comparison with the 
simple box. The evolution is a two-step process, quickly 
decohering kinematically to what would be the final level 
for the simple box, and then decaying exponentially to 
the final level relevant when the spectrum has a wing. 

For a spectrum with a long tail such as the Gaussian 
spectrum, the resonance always falls into the spectral 
range. In the sin 20 = limit one can not obtain sta- 
ble synchronized oscillations — eventually the system de- 
coheres, no matter how large /x, although the exponential 
decay is very slow for large fi. For nonvanishing mixing 
there exists a threshold so that for fi > fj,o partial syn- 
chronization arises. In Fig. [13] we show the Gaussian 
RAifJ-) for several mixing angles. In contrast to the box 
spectrum, there is always a threshold and it moves to 
larger fi as expected. 




Time 

FIG. 12: Evolution of order parameter for box spectrum with 
sin 20 — 0.04. Top curve: Simple box {b = 0). Bottom curve: 
b = 0.01. Dashed lines show the asymptotic values. 



1.0 




FIG. 13: Asymptotic order parameter for Gauss spectrum 
with mixing angles sin 26 = 1, 0.1 and 0.01 (left to right). 



VII. CONCLUSIONS 

We have investigated the most elementary case of col- 
lective neutrino oscillations, an isotropic gas prepared in 
a single flavor with a fixed density. The seemingly simple 
question of partial synchronization for a large but flnite 
density reveals a rich variety of possible solutions. 

The idea of a late-time asymptotic solution has led us 
to formulate self-consistency relations that can be evalu- 
ated in the sudden approximation and that provide the 
asymptotic solution from the spectrum alone. We flnd 
excellent agreement of the asymptotic solutions such pre- 
dicted with those from a numerical solution of the full 
EoMs. The agreement persists even when the true solu- 
tion does not resemble the sudden approximation and we 
suspect that these relations are exact. 

Partial synchronization requires a minimum density (a 
minimum effective interaction strength fi). For smaller 
densities, decoherence is complete. For larger /x, syn- 
chronization is partial and becomes perfect for /i — > cxi. 
This "phase transition behavior" derives from a simple 
model where one coarse-grains the spectrum in terms of 
two polarization vectors, representing two halves of the 
spectrum. The system is represented by a particle mov- 
ing in a double-well potential, and depending on /i it is 
trapped in one well or traverses both. 

Non-maximal mixing enables a parametric resonance 
where those polarization vectors near the resonance fre- 
quency are driven away from their initial orientation. 
The back-reaction on the system consists of decoherence. 
This effect can only happen if the spectrum actually in- 
cludes the resonance frequency, but this is always the case 
for spectra with long tails. The exact outcome depends 
on the spectrum and the mixing angle. 

Decoherence can thus arise in different ways. In the 
dilute limit it consists of purely kinematical de-phasing 
of different modes. In the interacting system the effect 
is similar, except that each mode precesses around a dif- 
ferent direction and this direction changes as the system 
decoheres. The self-induced resonant decoherence effect 
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that we have discovered, on the other hand, is dynamical 
in the sense that an internal resonance of the system leads 
to an exponential decay of the global polarization vector. 
This decay can be slow, depending on parameters, and 
introduces a new time scale. 

Our study is complementary to Pantaleone Q who 
used the completely incoherent state as an initial condi- 
tion. Based on an analytic argument, Pantaleone found 
that such a system will never develop spontaneous po- 
larization ("neutral stability"). He also argued that a 
collective mode can be supported for any strength of /i, 
in seeming contrast to our finding of a minimum required 
fi to maintain partial coherence. However, this ques- 
tion depends on the assumed initial condition. Beginning 
with all polarization vectors aligned with each other (our 
case), indeed a minimum /i is required for a coherent 
motion to survive. On the other hand, one can set up a 
pure precession mode for any strength of fi where all po- 



larization vectors are arraiiged in a co-rotating plane, not 
along a common direction [isl [iGj . In practice this can be 
achieved by beginning with — >■ oo (synchronized limit) 
and then reduce /i adiabatically to any desired value. 
This procedure leads to a spectral split if fi is reduced all 
the way to zero, and to a pure precession mode for any 
non- vanishing fi. 
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